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ABSTRACT 

Recent work has suggested that the fraction of obscured AGN declines with increasing luminosity, 
but it has been difficult to quantify this trend. Here, we attempt to measure this fraction as a 
function of luminosity by studying the ratio of mid-infrared to intrinsic nuclear bolometric luminosity 
in unobscured AGN. Because the mid-infrared is created by dust reprocessing of shorter wavelength 
nuclear light, this ratio is a diagnostic of /obsc, the fraction of solid angle around the nucleus covered 
by obscuring matter. In order to eliminate possible redshift-dependences while also achieving a large 
dynamic range in luminosity, we have collected archival 24 micron MIPS photometry from objects with 
z~l in the Sloan Digital Sky Survey (SDSS), the Great Observatories Origins Deep Survey (GOODS) 
and the Cosmic Evolution Survey (COSMOS). To measure the bolometric luminosity for each object, 
we used archival optical data supplemented by GALEX data. We find that the mean ratio of 24 fj,m 
to bolometric luminosity decreases by a factor of ~3 in the Lb o i=10 44 -3xl0 47 ergs s" 1 range, but 
there is also a large scatter at constant L\, \. Using radiation transfer solutions for model geometries, 
we show how the IR/bolometric ratio relates to /obsc and compare these values with those obtained 
obtained from samples of X-ray selected AGN. Although we find approximate agreement, our method 
indicates somewhat higher values of /obsc, particularly in the middle range of luminosities, suggesting 
that there may be a significant number of heavily obscured AGN missed by X-ray surveys. 
Subject headings: galaxies: active - quasars: general - infrared: galaxies 



1. UNIFIED AGN SCHEMES AT LOW LUMINOSITY 



Foll owing the initial recognition by lAntonucci fc M iller 
(1985) that the nucleus of the prototypical type 2 Seyfert 
galaxy NGC 1068 must be surrounded by a toroidal 
belt of gas and dust, a tremendous amount of evi- 
dence has accumulated in support of the idea that sim- 
ilar obscuring tori surround the nuclei of ot her type 2 
Seyfe r t galax i es (as reviewed, for example, in lAntonuccil 
119931: iKrolikl I1999D . Spectropolarimetry (the method 
originally employed by Antonucci and Miller) applied 
to both type 2 Seyfert galax ies and radio galaxies 
(|di Serego Alighieri et al.l Il994h often reveals in polar- 
ized light broad optical/UV emission lines and a strong 
non-stellar continuum that are essentially invisible in the 
total flux spectrum. In many examples of both these 
types of AGN, conical regions of bright line emission 
from highly-ion ized elements point at the otherwise ob - 
scured nucleus (|Ferruit et alJ 120001 : ISchmitt et al.ll2003f ). 
Similarly, both type 2 Seyfert galaxies and radio galax- 
ies often exhibit X-ray continua that either show evi- 
dence for very large column densities of material along 
the line of sight or are extremely weak, l ikely indicat- 
ing that the obscuration i s Compton thick (|Risaliti et al.l 
U999tlTreister et al.ll200"l . In two cases, NGC 1068 and 
the Circinus galaxy, the obscuring torus c an be seen 
direc t ly in interferometric infrared imaging (|Jaffe et al.l 
120041: iTristram et al.|[2007t) . 

The principal effect of this toroidal obscuration is that 
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observers trying to see the nucleus along lines of sight 
that pass through it are prevented from seeing any- 
thing but dust-reprocessed infrared continuum and (per- 
haps) hard X-rays. Consequently, most of the classic 
signatures of AGN — broad optical/UV emission lines, 
strong optical/UV non-stellar continuum, strong X-ray 
continuum — are obliterated when the object is seen from 
such a direction. Only when the line of sight passes 
through the central hole of the torus can all these fea- 
tures be seen, and the object is perceived as a type 1 
AGN. 

In the nearby Universe, where we can see large numbers 
of low luminosity AGN, large statistical samples have 
been assembled in order to measure the luminosity func- 
tions of both type 1 and type 2 AGN. Ideally compared at 
matched bolometric luminosity, but more often in terms 
of their luminosity in a single band or feature, the ratio 
of their numbers can be immediately interpreted as the 
ratio of unobscur ed to obscured solid a ngle. For example, 
lHao et ail (|2005h and ISimpsonl (|2005h . using samples of 
optically-selected AGN from the Sloan Digital Sky Sur- 
vey, found that the ratio of type 1 objects to type 2 
objects increases with increasing luminosity, from ~ 2/3 
at an [OIII] 5007 luminosity of 10 6 L Q to ~ 2 for line lu- 
minosi ties ~3x 10 7 L^. U sing s amples of X-r a y sele cted 
AGN, iSteffen et all (pool and iBarger et all ((2005h ar- 
gued for a similar shift in the ratio of unobscured to 
obscured, suggesting that the unobscured variety come 
to dominate the total population for 2-8 keV luminos- 
ity <~, 10 44 erg s" 1 . Similar result s were obtained by 
lUeda et al.l (|2003f ). lLa Franca et al.l (|2005l ) and others, 
also using X-ray selected AGN samples. In fact, hints of 
such a trend were already see n in the much smaller 3CR 
radio sample iLawrencel ()1991f ) . 

However, there remain significant uncertainties in the 
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measured population ratio at high luminosity. For exam- 
ple, Compton-thick AGN are completely missing from X- 
ray-selected samples. Other sample construction meth- 
ods (e.g., the very large sample of type 2 quasars se- 
lec ted from the SPSS ba sed primarily on [OIII] emission 
bv lZakamska et aT1l2003f l can potentially yield quantita- 
tive estimates of comparative luminosity functions, but 
require substantial work in order to turn catalogs into 
space densities. Although much effort has been made 
in the last few years to construct AGN samples from 
mid-infrared observations with the Spitzer Space Tele- 
scope, the selection effects remain too poorly understood 
and the sample size is sometimes too small to permit 
extraction of a trend in the type ratio as a function of 
luminosity (IStern et al.ll2005tlAlonso-Herrero et al"l2006t 
iMartmez-Sansigre et alj|2006t lLacv et al.ll2007f ). 

Here we attempt a different approach to the prob- 
lem of measuring the unobscured/obscured ratio: us- 
ing the ratio of reradiated mid-infrared continuum to 
bolometric luminosity in type 1 AGN a s a surrogate 
(cf. the suggestions in iLawrencel 119911 : iBarger et al.l 
2005). Because the obscuration transforms essen- 
tially all the nuclear radiation incident upon it into 
mid-infrared continuum via dust reradiation, and de- 
tailed radiation transfer calculations suggest that most 
of that continuum is radiated toward observers in 
the type 1 direction no matter what the d e tailed - 
arrangement of dust may be dPier fc Krolikl 119921: 
Granato & Danese 1994; Efstathiou & Rowan-Robinson 
1995 1 : iNenkova et al.l 120021: Ivan Bemmel fc Dullemondl 
20031 iDullemond fc van Bemmelll2005l ). this ratio should 
give a good indicator of the ratio of unobscured/obscured 
solid angle. The relative ease of identifying type 1 AGN 
across a wide dynamic range in luminosity also gives this 
method a number of advantages relative to those depen- 
dent o n actually counting type 2 objects. iMaiolino et al.l 
(2007) have recently made an effort to implement this 
program, but the details of their method differ in signif- 
icant ways from ours; we will comment on specific con- 
trasts as they arise. 

This paper is structured as follows: In section §2 
we present the selection criteria and basic observational 
properties of the sources used in this work. In section §3 
we study the correlations found in this sample, while the 
significance and interpretation of these trends are dis- 
cussed in §4. Our conclusions are presented in §5. When 
required, we assume a ACDM cosmology with /io=0.71, 
^™=0.3 and £1^=0.7, in a greement with the m ost recent 
cosmological observations (jSpergel et alj |2007) . 

2. SAMPLE DEFINITION AND DATA SOURCES 

In order to eliminate any possibility that what we find 
involves an evolutionary effect rather than a luminos- 
ity dependence, we chose a sample of AGN that all have 
nea rly the same redshift ( 0.8 < z < 1.2), in contrast with 
the IMaiolino et alj (|2007f) sample which includes sources 
at all redshifts. This particular redshift is a convenient 
choice because it is large enough that luminous quasars 
were common, but not so large that low luminosity AGN 
are too faint to detect. To achieve the greatest possible 
dynamic range in luminosity, we made use of three sam- 
ples with complementary properties: one that has very 
large solid angle but is relatively shallow (Sloan Digital 
Sky Survey: SDSS), one that has smaller solid angle but 



goes deeper (Great Observatories Origins Deep Survey: 
GOODS) and one intermediate approach (Cosmic Evo- 
lution Survey: COSMOS). 

For our infrared band, we used the Spitzer Space Tele- 
scope MIPS 24/zm channel. The corresponding rest- 
frame wavelengths are ~12/mi; by requiring z<1.2, we 
avoid most of the contamination from the 9.7/xm sili- 
cate feature. For the highest redshift sources, this sil- 
icate feature can influence our 24//m photometry, how- 
ever this will affect only a small fraction of our sample. 
In addition, at this relatively short wavelength, stellar- 
heated dust is usually relatively faint. From the compos- 
ite QSO spectrum with individually detected PAH fea- 
tures of Schweitzer et al. (2006; Figure 2 top panel), we 
can estimate the fraction of integrated flux contributed 
by the 11.3 /im PAH feature in the MIPS 24 micron band 
at 0.8<z<1.2, and we find a maximum contribution of 
only 1.3%; hence this contribution is completely negli- 
gible. Moreover, because we use only unobscured AGN 
with Lboi from 10 44 to 10 47 5 ergs s _1 , stellar heating is 
unlikely to be a significant contaminant. The AGN from 
the SDSS were selected based on their optical properties, 
while in the case of the GOODS and COSMOS sources 
they were X-ray selected. In any case, the selection is 
completely independent of their infrared properties, so 
it will not bias our results. In only one case, a SDSS 
source is located in the field of view of a Spitzer observa- 
tion, but a counterpart was not detected. In that case, 
we estimate the corresponding upper limit as described 
below. 

For our indicator of the nuclear luminosity, we used 
optical and UV (GALEX) data to estimate the bolomet- 
ric luminosity individually for each source. Because the 
GALEX NUV band (1750-2800 A ) corresponds to 875- 
1900 A in the rest-frame, we are able to measure directly 
the flux in much of the spectral peak region. The details 
of our method are presented in §2.5. In addition, we have 
also found it is sometimes convenient to use the SDSS i- 
band luminosity as a standard of comparison because in 
our redshift range it corresponds to rest-frame B-band. 

2.1. SDSS Sample Data Collection 

In order to construct a large sample of high luminosity 
unobscured quasars, we us ed the results fro m the Sloan 
Digital Sky Survey (S DSSJYork et al]l200Ch. spe cifically 
from data release 5 ()Adelman-McCarthv et ail I2006) 4 . 
Only sources classified as quasars by the SDSS pipeline, 
and thus presenting high- ionization broad emission lines, 
were considered. A detailed description of the SPSS 
quasa r classification scheme was presented bv lYip et al.1 
(|200l . There are 11937 SDSS/DR5 quasars in the 
0.8<£<1.2 redshift range. 

Spitzer data for these SDSS quasars were obtained 
from the archive using the Leopard v6.1 software 5 . A 
total of 638 Astronomical Observation Requests (AORs) 
were found with at least one source in our sample in 
the field. In many cases, these correspond to different 
observations of the same field, so only 206 quasars in 
our SDSS/DR5 sample were observed by Spitzer/MIPS. 
Only the post-BCD (Basic Calibrated Data) frames for 

4 Data available at http: / /ww w. sdss.org/dr5] 

5 Leopard can be obtained at 

http: / / ssc. spitzer . caltech.edu /propkit /spot/ 
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each AOR were downloaded. The Mosaicking and Point - 
source Extraction (MOPEX; iMakovoz fc Marleaull2005t ) 
package version 030106 6 was used in order to extract the 
sources and calcul ate fluxes from each APR. The proce- 
dure described bv lMakovoz &: Marlea u (2005) 7 , was fol- 
lowed. Briefly, we did a first pass extraction including 
the brightest sources only, in order to calculate the point 
response function (PRF) for each mosaic. Using this 
specific PRF we then extracted sources with a signal- 
to-noise higher than 5. In the last step, the flux for each 
source was estimated by performing PRF fitting. In or- 
der to check this procedure, we downloaded the data from 
the Great Observator ies Origins Deep Survey (GOODS; 
iDickinson et al.ll2003f ) and compared the results obtained 
by following this procedure with those reported by the 
legacy survey team 8 . The flux densities presented in the 
GOODS catalog were converted into fluxes assuming a 
MIPS-24 bandwidth of 4.7 /um, as reported by the MIPS 
data handbook. We found very good agreement between 
the two independent studies. On average, the fluxes re- 
ported by the GOODS team were higher by ~ 10% with 
a spread of less than ~10%, even after incorporating an 
aperture correction. Given that this small offset will not 
affect our results, we do not apply any correction to our 
fluxes. 

The SDSS sample was then matched to the sources 
found in the Spitzer archive, using a maximum search 
radius of 10". This maximum nominal separation was 
chosen both to minimize the number of false counterparts 
and to take into account the positional uncertainties of 
both Spitzer and SDSS. Altogether, we found in these 
Spitzer images counterparts for 205 sources at a signal 
to noise higher than 5. In the remaining case, SDSS 
J024256. 93-001558.0, a very faint counterpart is visible 
on the Spitzer image, but not significantly detected so we 
calculate a 3-er upper limit from the background standard 
deviation measured around the SDSS position. The dis- 
tribution of distances between the SDSS position and the 
Spitzer counterpart can be seen in Figure Q] The average 
separation between the optical position and the IR coun- 
terpart was ~0.5". For each SDSS source, we calculated 
the probability of a false match with a 24 /im counter- 
part, assuming a Poissonian distribution of sources and 
the source density of the corresponding AOR. We found 
that the maximum probability of a false match for a given 
source was ~ 1.3%, with an average probability of 0.02%. 

The observed properties of these 206 sources can be 
found in Table [U These objects were on average our 
most luminous, with bolometric luminosities ranging 
from 2 x 10 45 erg s" 1 to 3 x 10 47 erg s _1 . Optical mag- 
nitudes in the AB system were obtained directly from 
the reported SDSS values. In all cases, the PSF fitting 
magnitude was used, as suggested by the SDSS team for 
sources that are unresolved at ground-based resolution. 
Although the errors in the SDSS photometry are negli- 
gible, since these are all relatively bright optical sources, 
errors in the MIPS-24 /xm photometry are estimated to 
be ~10-20%, including both measurement and system- 

6 MOPEX can be downloaded from 

http: / /ssc. spitzer. caltech.edu/postbcd/download- mopex.html 

' A more detailed description can be fo und at 
http: / / ssc. spitzer . caltech.edu / postbcd / running- prf . html 

A The GOODS Spitzer data can be found at 
http: / / data.spitzer.caltech.edu/popular / goods/ 
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Fig. 1. — Distribution of distances between the optical posi- 
tion from the SDSS and the Spitzer MIPS counterpart. Average 
separation is 0.5". 



atic (~ 5% according to the MIPS data handbook) un- 
certainties. 

2.2. GOODS Sample Data Collection 

In order to study lower luminosity sources, we include 
sources detected in both the north and south GOODS 
fields. A summa ry of the properties of these sources 
was presented by iTreister et al.l (|2006f) . Our selection 
criteria were the same as for the SDSS sample: unob- 
scured sources and 0.8<z<1.2. In this case, a source 
was classified as unobscured if broad lines were present 
in the optical spectrum and the equivalent hydrogen col- 
umn density of obscuration Nh in the X-ray spectrum 
was smaller than 10 22 cm" 2 . With these criteria, we 
found a total of 10 sources, 4 from the GOODS-N and 6 
from the GOODS-S fields. These sources are in gen- 
eral intrinsically fainter than the SDSS sources, with 
bolometric luminosities ranging from lxlO 44 erg s _1 to 
1 x 10 46 erg s" 1 . 

Optical magnitudes for these sources obtained from the 
Hubble Space telescope with the Advanced Camera for 
Surve ys (ACS) imager were presented by the GOODS 
team (Giavali sco et al.ll2004l ) 9 . In order to avoid prob- 
lems related to the different spatial resolution of the 
GOODS and the SDSS images, we used aperture magni- 
tudes with a 2" diameter to calculate the optical fluxes 
of the GOODS sources. 

The Spitzer properties and redshifts of these sources 
were obtained from the compilation of ITreister et al.l 
(pOOl . and are summarized in Table [2] The ID num- 
bers in t able [H correspond to t he X-ray identifications 
from the lAlexander et al.l (|2003f ) catalog. Flux densities 



' Available at http://archive.stsci.edu/pub/hlsp/goods/catalog_rl/ 



4 



Treister, Krolik and Dullcmond 



in the 24 /zm band were converted into fluxes by multi- 
plying them by the bandwidth of 4.7 //m. 

2.3. COSMOS Sample Data Collection 

The COSMOS survey, covering a total area of 2 deg 2 
at higher flux limits than GOODS, provides sources 
that fill the gap in luminosity between the SDSS and 
GOODS sources. Their span of luminosities is from 
2 x fO 44 erg s" 1 to 1 x 10 46 erg s~ 4 . We selected 
sources from this survey based on the optical spec- 
troscopy observations of the XMM X-ray sources us- 
ing the M agellan/IMACS mult i-object spectrograph pre- 
sented by iTrump et al.1 (|2006l ). Using our selection cri- 
teria (0.8<z<1.2 and broad optical emission lines) we 
obtained a sample of 19 sources. 

The COSMOS field was completely observed by Spitzer 
using both IRA C and MIPS as part of the S-COSMOS 
legacy program (|Sanders et al. 2007) . Using the public 
catalog of Spitzer sources generated by the COSMOS 
team 10 , we detected 14 of the 19 X-ray selected sources in 
the MIPS 24-/im band. The observed properties of these 
sources are presented in Table |3l As for the GOODS 
sources, flux densities in the 24 /jm band were converted 
into fluxes by multiplying them by the bandwidth. 

In summary, our sample consists of 230 sources span- 
ning a range in bolometric luminosity from Lboi=10 44 to 
lO 47 5 erg s^ 1 . All the sources in our sample have mea- 
sured spectroscopic redshifts and present broad emission 
lines. In the IR, all these sources are bright, with lumi- 
nosities from Ltr=10 43 5 to 10 46 ergs s _1 . 

2.4. G 'ALEX Data 

As a large fraction of the bolometric luminosity in un- 
obscured AGN, ~15-20%, can be found in the UV (e.g., 
Elvis et all [l994t iRichards et ail l2006t iTrammell etafl 
20071) . it is worthwhile to look for UV counterparts of our 
sources in data taken by the Galaxy Evolution Explorer 
(GALEX: [Martin et al.ll2005D . Specifically, we made use 
of data release 3 11 . In the case of the SDSS sources, 
we looked for UV counterparts in a 10" radius around 
the optical position. The median separation over the 
whole sample was much less than the maximum permit- 
ted, ~0.5", and the actual maximum separation was only 
3.5". 170 of the 206 sources in the SDSS sample were de- 
tected in at least one of the GALEX bands, for a detec- 
tion rate of 83%. Of these 170 sources, 61 were detected 
only in the reddest NUV band and 4 detected only in the 
FUV band. The GALEX properties of the SDSS sources 
are presented in Table [TJ 

GALEX coverage of both GOODS fields is signifi- 
cantly deeper. Whereas the all-sky survey exposure time 
was only 0.1 ksec, the North field is being observed for 
100 ksec as part of the deep imaging survey, and the 
South field is being observed for 200 ksec as part of 
the ultra-deep imaging survey. At the time of data re- 
lease 2, the total time spent on the South field was only 
~76 ksec, while ^95 ksec had been accumulated in the 
North field. Consequently, only 4 of the 6 sources in 
the South field were detected by GALEX, but all the 
sources in the North field have a UV counterpart. In Ta- 
ble El we present the GALEX properties of the GOODS 



10 Catalog available at http:/ /irsa. ipac.caltech.edu/Missions/cosmos 

11 Available at http://galex.stsci.edu/GR2/ 



sources. The entire COSMOS field is one of the tar- 
gets of the GALEX Ultra-Deep Imaging Survey. Hence, 
it is not surprising that all the sources in our sample 
with Spitzer-MIPS detections were also detected in the 
GALEX observations. The UV properties of the COS- 
MOS sources are presented in Table [3] 

2.5. Bolometric Luminosities 

Because we believe the infrared continuum to be the 
result of reprocessing shorter wavelength radiation from 
the nucleus that strikes the dusty torus, including it in 
our estimate of the intrinsic nuclear bolometric luminos- 
ity when we have an unobscured view of the nucleus 
would be double-counting. Instead, we do the best we 
can to sum all flux from wavelengths shorter than ~l/xm; 
the exact long wavelength cut-off is not important be- 
cause the integral is in general dominated by the UV. 
We neglect X-ray c ontributions, but the se are generally 
small, < 10% (e.g., IRichards et al.ll2006t) . 

In order to estimate the bolometric luminosities, we 
combined the photometric information from GALEX 
NUV (1750-2800A ) to z-band (~l/xm). In the opti- 
cal range, where there are only very small gaps in wave- 
length between the available photometric bands, we ap- 
proximate the integrated luminosity by the sum of the 
observed luminosity in each band. For UV fluxes, we em- 
ploy GALEX data wherever possible. Specifically, we use 
the GALEX NUV channel, which for our sample trans- 
lates to rest-frame wavelengths 875-1400 A . Whenever 
the object was either not observed by GALEX, or has 
only an upper bound, we estimate the flux in the NUV 
band by extrapolating from the bluest filter available (it- 
band for COSMOS and E-band for GOODS) using the 
median quasar SED from lRichards et al.l ((2006j).Botli be- 
cause GALEX data are missing from only a small frac- 
tion of our sample and because th e intrinsic dispersion i n 
bolometric corrections is ~ 50% ([Richards et al.l 12006). 
this extrapolation should not bias significantly our re- 
sults. To account for the gap between the u band (rest- 
frame ~ 1750 A ) and the GALEX NUV band, we lin- 
early interpolate in F v . At even shorter wavelengths, 
which are often not accessible in direct observations, 
we assume a power-law spec trum (f v oc y a ) wi th slope 
a = -1.76, as reported bv iTelfer et~aT1 (|2002f ) in the 
500-1200 A regime. The amplitude of this power law 
is fixed by matching to either the GALEX photometry 
(when available) or the extrapolated flux from the SDSS 
composite spectrum when G ALEX data are not a vail- 
able. From the templates of IRichards et all (|2006f ). we 
estimate that ^20% of the bolometric luminosity is emit- 
ted in this wavelength range. 

The approach described above incorporates all the 
available information in the o ptical-UV range. In 
contrast, iMaiolino et alj |2007) computed bolometric 
luminosities by multiplying the observed rest-frame 
5100 A luminosity by a fixed bolometric correction. 

3. TRENDS 

Our central result is shown in Figure El There we plot 
the ratio of MIPS 24 /im luminosity, LmipSj to our esti- 
mate of the bolometric luminosity, Lboi, of each individ- 
ual object. Two properties of the data are immediately 
htajjparent from this figure: that the mean value of this ra- 
tio decreases by order unity over the factor of 300 range 
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Fig. 2. — Ratio of 24 fim to bolometric luminosity as a function 
of bolometric luminosity for all the sources in the sample. The 
open circles show the ratio for sources not detected by GALEX. 
Triangles with error bars show the average value in each luminosity 
bin, with 1-a vertical error bars. 



in luminosity these data span; and that at any single 
luminosity there is a very large dispersion. 

Although (iMiPs/i'boi) does change over this range of 
bolometric luminosities, it is not by a large factor: a 
factor of 3 is the best estimate. This figure may be com- 
pared wit h the factor of 2 contr ast given by the fitting 
formula of Maiolino et al. (2007), meant to describe the 
sample means from ALa(5100 A ) = 10 43 ergs s _1 to 
10 47 ' 5 ergs s _1 , a somewhat greater dynamic range than 
spanned by our sample. Nonetheless, our sample size 
is large enough for the bin-means to be quite narrowly- 
defined, and the trend is clearly significant. Several sta- 
tistical tests support this conclusion. To evaluate the sig- 
nificance in difference between the mean values in the two 
bins with the largest number of sample points (the sec- 
ond and third highest in luminosity), we apply Student's 
t-test; it shows that they are different at the 0.1% signifi- 
cance level, while the means of the highest and the third 
highest luminosity bins are different at the 10 -5 level. 
Similarly, the Spearman rank-correlation test finds that 
the correlation between log(Lboi) and log(LMiPs/-Lboi) is 
significant at the 10 -7 level. A conventional Pearson cor- 
relation test fails on these data because there are so many 
more SDSS quasars than COSMOS and GOODS sources, 
and the contrast in the mean infrared/bolometric lu- 
minosity across the narrower range of luminosity where 
most of the SDSS quasars are found is relatively small. 

At the same time, the dispersion is very large. Consid- 
ered in logarithmic terms, the rms fluctuation within the 
bins is 0.25-0.35, equivalent to a multiplicative factor of 
1.8-2.2, comparable to the contrast in the mean across 
our entire luminosity range. 
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Fig. 3.— Same as Figure [2] but showing the GALEX NUV to 
bolometric luminosity (left panel) and i-band to bolometric lumi- 
nosity (right panel) ratios as a function of bolometric luminosity. 
Dashed error bars show the mean values in the lowest luminosity 
bin if source COSMOS J100129. 83+023239.0 is excluded. 



That this trend (and large dispersion) is not related to 
effects in the intrinsic spectra of AGN can be seen by con- 
sidering two other flux ratios as functions of Lboi (Fig-EJ)- 
In the mean, there is no significant change in the ratio of 
either the i-band luminosity or the NUV-band luminosity 
to bolometric. The small change observed in the z-band 
to bolometric luminosity ratio is heavily influenced by 
the single source COSMOS J100129.83+023239.0, which 
has a ratio ^8x higher than the rest of the sources in 
that bin. If this source is removed, the average in the 
lower luminosity bin decreases from 0.052 to 0.041 and 
the trend disappears, as can be seen in Fig. [3] Performing 
a Spearman test, we found that the L^/Lboi correlation 
with Lboi is ^3 orders of magnitude less significant than 
the I/MiPs/iboi correlation and has roughly the same sig- 
nificance as the ^Nuv/^boi correlation. 

Moreover, there is far less intrinsic scatter at fixed -Lboi 
in either of these ratios than in the infrared to bolomet- 
ric ratio. The standard deviation of LiR/Lboi in each 
luminosity bin is a factor of ~2.3, while for Lj/Lbolj it is 
~1.5 and for -LNuv/-Lboi it is 1.12. The magnitude of the 
error induced by using a fixed bolometric correction to es- 
timate bolometric luminosities from a given optical band 
is shown by the dispersion of the points plotted in both 
panels of Figure [31 Use of the i-band and a fixed bolo- 
metric correction, for example, entails an uncertainty in 
the bolometric luminosity of ~ 50%. 

4. SIGNIFICANCE 

4.1. Evolution? 

Just as it is possible that the physics controlling the 
geometry of obscuration is correlated in some way with 
the bolometric luminosity of the AGN, so too, it may 
change as a function of the age of the Universe. There is 
evidence, for example, that at fi xed luminosity, .fqbsc in - 
creases with increasing redshift (jTreister fc Urrvl [2006*1: 
this suggestion, however, remains controversial (e.g., 
iGilli et aLll2007| ). Factoring out any possible redshift- 
dependence was a strong motivating factor for our choice 
of a sample in which all the objects have approximately 
the same z. 

In fact, the precise character of the correlation between 
/obsc and Lboi also enters into the debate over the pos- 
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sible dependence of /obsc on z. Currently-available sam- 
ples that probe to high redshift are in general strongly 
incomplete at low luminosity because it is so difficult to 
obtain spectra for objects that faint. If one's estimate of 
the obscured fraction at high redshift is based solely on 
the detectable sources (which are, of course, high lumi- 
nosity objects), it will be biased to a low value if / b sc 
systematically declines with increasing luminosity at all 
redshifts. 

4.2. Relation between iMiPs/^boi and the obscured 
solid angle 

As we explained earlier, the ratio of infrared repro- 
cessed light to bolometric luminosity should increase as 
the fraction of solid angle occupied by the dusty obscur- 
ing gas /obsc increases. If the reprocessed IR were radi- 
ated isotropically, the ratio Lmips / -^boi would be propor- 
tional to /obsc- It is much more likely, however, that the 
reradiation is anisotropic. Consequently, the relationship 
between these two quantities is not necessarily precisely 
linear, nor does it directly yield /obsc wit hout the intro- 
ducti on of any additional parameters (cf. Maiolino et al. 
2007). At the very least (i.e., in the event of isotropic 
radiation), there is a proportionality constant to be de- 
termined. 

At present, there are two lines of evidence regarding 
the directionality of the reradiated IR. On the observa- 
tional side, the re are suggestions that it may not b e far 
from isotropic (jLutz et ail 12004 iHorst et alj|2007l ). In 
both of these studies, it is claimed on the basis of small 
samples of AGN that the mean ratio of mid-IR flux to 
intrinsic hard X-ray flux is the same for both obscured 
and unobscured objects. However, in both of these stud- 
ies, the intrinsic hard X-ray luminosities of the obscured 
objects are in general substantially smaller (by at least 
an order of magnitude) than those of the unobscured 
objects. Given the correlation we have found between 
^Mips/^boi and iboi, one might then have expected that 
in their samples the mid-IR to X-ray flux ratio would be 
even higher for the obscured than for the unobscured if 
the IR radiation were isotropic. We therefore do not re- 
gard these studies as making a strong case for isotropic 
radiation. 

On the other hand, there are strong physical and 
theoretical arguments for thinking that the radiation 
is anisotropic. The obscuration is, after all, strongly 
aspherical — if it were otherwise, there would be very 
few unobscured objects to see. If the torus is opti- 
cally thick in the radial direction in the mid-IR and that 
optical depth is larger than the vertical optical depth, 
most of the mid-IR flux should escape roughly paral- 
lel to the torus axis, with rather less escaping in the 
equatorial direction. Every detailed radiation transfer 
model published shows this tendency, although the pre- 
cise contrast between the 12/im flux toward the pole 
and toward the equator varies depe nding on the details 
of the adopted density d i stribution dPier fc Kroliklll993l: 
Granato & Dancsc 1994; Efstathiou & Rowan-Robinson 



19951 iNenkova et al.l l2002t Ivan Bemmel fe Dullemondl 
200alDullemond fe van Bemmelll2005D . The most recent 



paper (which presents both smooth and clumpy models) 
finds ratios between the total infrared flux in the equa- 
torial plane and at 20° from the axis to vary between 
0.06 and 0.3, depending on the specific model. Making 



the approximation that very little flux is directed toward 
observers in the obscured solid angle, we expect that the 
MlPS-band luminosity should scale relative to the bolo- 
metric luminosity approximately as 



dLuiPS 

dn 



/l2(0) ; 



^bol 



/. 



obsc 



47T 1 - f ohB 



(1) 



Here /12(d) is the fraction of the total dust-reprocessed 
luminosity falling within the MIPS band (~ 11-13/im 
in the rest-frame at z ~ 1). Even within the range of 
viewing angles 9 that provide an unobscured view, we 
may expect order unity variations in this quantity. 

To confirm this scaling and estimate f 12(d), we stud- 
ied the results of a number of model calculations. These 
models were performed with t he code described in 
iDullemond fc van Bemmei (2005). For all of them, the 
obscuration was assumed to be confined to within a con- 
stant opening-angle spherical wedge whose outer radius 
was 30 times greater than its inner radius. We consid- 
ered two density profiles: a single constant value and a 
radial decay cxr -1 . For both density models, when the 
torus is optically thick in the mid-infrared, the approxi- 
mate model given by equation 1 is confirmed provided the 
opening angle is at least ~0.7 radians. Crudely speaking, 
the solid angle-weighted mean value of /12 (for viewing 
angles permitting an unobscured view of the nucleus) 
does not vary a great deal with torus opening angle: it 
varies from ~0.06 to ~0.08 depending on the model de- 
tails for the same range of opening angles for which equa- 
tion 1 holds. When the opening angle becomes smaller 
than this, /12 can begin to diminish, with the amount 
sensitive to whether the torus is truly a spherical wedge 
or has a flatter outer envelope. At fixed opening angle, 
/12 does vary somewhat with viewing angle, typically be- 
ing perhaps 50% greater than the mean for nearly polar 
viewing angle and dropping sharply when the viewing an- 
gle becomes close to the torus opening angle. The latter 
effect is also sensitive to the specific geometry (spheri- 
cal wedge) adopted in these models, which also depends 
slightly on the wavelength of the emission. 

The ratio of IR flux to bolometric seen by an observer 
in the unobscured direction is then 



£mips 
Lbo\ 



f 12(e) 



obsc 



1-/ 



obsc 



(2) 



To this level of approximation, the obscured fraction as- 
sociated with a given value of £iviiPs/£boi is 



/< 



obsc 



where 



l + q' 
f 12(0) 

LlAlVs/Lhol 



(3) 
(4) 



Using the solid angle-weighted means of /12 taken 
from our models, we can transform the bin-means of 
iMiPs/^boi shown in Figure [2] to estimates of /obsc- The 
results are shown in Figure [4] It appears that the mean 
value of /obsc falls from ~ 0.9 at the lowest bolometric 
luminosities in our sample to ~ 0.3-0.4 at the highest 
luminosities. We strongly caution, however, that there 
may be a systematic error of order unity in (/12} due to 
the specific assumptions made in those models regarding 
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Fig. 4. — Two different inferences of / bsc : black curve with error 
bars from X-ray data, as described in the text; triangles with error 
bars from the method of this paper assuming values of /i2=0.06 
(squares; red in electronic version) and 0.08 (triangles; blue in elec- 
tronic version; displaced by logLbol=0.1 for clarity). The horizon- 
tal error bars simply describe the width of the luminosity bin; the 
vertical red and blue error bars are dashed to indicate that they 
are model-dependent in the sense discussed in the text. Their ex- 
tent was determined by translating the error bars of Fig. [2] using 
Eqn.[3] 



the density distribution within the envelope that deter- 
mines /obsc- A variety of parameters (e.g., the ratio of 
outer radius to inner, the degree of inhomogeneity and 
clumping, etc.) may quantitatively alter (fi2)- 

We compare this dependence of the mean /obsc as 
a function of Lboi with the relation that has been 
found in X-ray studies in Fig ure ffl For thi s figur e 



we use the X-ray sample of iTreister fc Urrvl (|2006l ) 
with the X-ray luminosity transformed to Lboi us- 
ing the luminosity-d ependent bolometric corrections of 
iMarconi et al.l (|2004l ). The level of agreement between 
the two independent methods of inferring /obsc is sur- 
prisingly good, considering the systematic errors in 
both. Most notably, X-ray surveys completely miss ob- 
jects with column densities 3> 10 24 cm -2 , and there 
are indications that these objects may be num erous 
(|Martmez-Sansigre et all 120061 ; IPolletta et al.ll2006h . In- 
deed, if one wished to extend credence to the somewhat 
larger / bsc found by the infrared method at intermedi- 
ate luminosities, it might be explained by a somewhat 
larger proportion of Compton-thick objects within the 
population. It i s also interesting that, as pointed out by 
iMaiolino et al.l (|2007l ). the values of /obsc derived from 
X-ray surveys include the effects of obscuration by both 
dust and gas, while the IR-derived values include only 
obscuration by dust. 

To close this subsection, we remark on another dis- 
tinction that may be important in evaluating the mean- 



ing of the /obsc estimated from iMiPs/^boi- It is pos- 
sible for obscuring matter of smaller optical depth to 
intercept some of the nuclear light farther out in the 
host galaxy and reradiate it in the mid-infrared even 
though this matter is whol ly independent of the "ob- 
scuring torus" proper ( e.g., lKeellll980l : iMcLeod &: Riekd 
ll995tlRigbv et al.|[200g ). Molecular clouds and other gas 
concentrations in the host have at most modest column 
densities when measured in terms of X-ray absorption 
(generally < 10 22 cm -2 ), but such column densities can 
still absorb large portions of an AGN's optical/UV con- 
tinuum. Consequently, they qualify as contributing to 
/obsc both in the sense that they can block our view of 
the optical/UV light from the nucleus and, depending on 
their temperature, in terms of their contribution to the 
mid-infrared luminosity of the system. 

4.3. The dispersion at fixed luminosity 

One of the virtues of using SDSS quasars is that the 
sample size, particularly for 5 x 10 45 erg s _1 < -Lboi < 
1 x 10 47 erg s _1 , is large enough to be truly statistical. 
As a result, we can obtain a reasonably reliable estimate 
of the sample dispersion as well as its mean. As we have 
previously remarked, the dispersion is quite large. 

We see no measurement error that could contribute 
in a significant way to the dispersion in iMiPs/^boi at 
fixed Lboi- The errors in the photometry we use to com- 
pute the bolometric luminosity are small, ^2%; the un- 
certainty in the MIPS fluxes is ~10-20%, as detailed in 

§ 2 - L . 

Intrinsic variability could also contribute, but in a sim- 
ilarly minor fashion except perhaps in a very small num- 
ber of objects. The time interval between the GALEX 
measurements, the SDSS photometry, and the Spitzer 
measurements was gen erally a few years . O n that 
timescale, SDSS data (jVanden Berk et all I2004D show 
that quasars typically vary by 0.1-0.2 mag in the op- 
tical, with the larger variability usually occurring at 
shorter wavelengths and in lower luminosity objects. A 
small fraction vary by as much as ~ 0.5 mag on these 
timescales, still much smaller than our observed disper- 
sion relative to the mean. If the infrared is, as we assume, 
the result of thermal reprocessing by dust, it cannot vary 
substantially on timescales shorter than several to ten 
years. 



According to the results oflSchwcitz er et aTl (|2006f ). the 
contribution of star formation to the total light in quasars 
is ~ 30% at far-IR wave lengths, in agreement with the 
results of lShi et all (|2007l ). who found a contribution of ~ 
25% at 70 /xm and ~ 10% at 24 /xm. Hence, extrapolating 
linearly to the MIPS band, we expect a contribution of 
star formation to the total IR light of ~ 5% at rest- 
frame 12 /xm; we therefore do not expect our results to be 
significantly affected by star formation processes in the 
host galaxy. In addition, because fobsc is larger for lower 
luminosity sources than for higher luminosity sources and 
the nuclear infrared luminosity is fobscLboi, we do not 
expect the contribution from star formation to increase 
significantly even at lower luminosities. 

There must, therefore, be a genuine large dispersion in 
^Mips/iboi- One possible source of this dispersion is the 
variation of /12 with viewing angle. Again relying on the 
models described above, a total range of a factor of 2 3 
may be expected from this source alone. Compared with 
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our rms logarithmic dispersion of a factor of ~2, or the 
factor of 100 total range seen in the most-populated bins, 
the expected variations in viewing angle can account for 
a significant part of this dispersion, but likely not all of 
it. 

This last statement must, of course, be qualified by the 
condition that our models accurately estimate the range 
in J '12(0)- It is entirely possible that changes in the de- 
tailed density distribution, even while keeping / bsc fixed, 
might introduce further variations in f 12(G)- Although 
the expected large optical depth in the 12/xm band is 
likely to limit the sensitivity of f 12(d) to changes in the 
radial density profile or density inhomogeneities, there 
could certainly be small, but detectable, quantitative ef- 
fects. 

Finally, there may, of course, be additional dispersion 
due to variations in / bsc at fixed luminosity. It is entirely 
possible for parameters other than Lboi to influence the 
obscuration's solid angle. Among the obvious candidates 
are the mass of the central mass black hole and the stellar 
mass contained within the region of the obscuring matter. 
The varying impact of magnetic forces, clump collisions, 
or other mechanisms could also conceivably play a part. 

4.4. Implications for surveys 

Particularly since the launch of the Spitzer Space 
Telescope, it has become feasible — and worthwhile — to 
construct s amples of A GN b as ed on mid -in frared 



selection (lLacv et al 



iMart mcz-Sansigre et al 



2004| iStern et all [20051 : 
120061) . Our results here 
demonstrate that the translation between mid-infrared 
flux and bolometric flux has an intrinsic scatter of a 
factor of ~ 2 up or down, even for unobscured AGN. 
Because, as we have discussed, -LiR./Lboi is generally 
smaller in obscured AGN than in unobscured, the 
bolometric correction based on infrared flux for them 
should be an even larger factor, with perhaps greater 



scatter. 



5. CONCLUSIONS 



We presented in this paper an alternative way to esti- 
mate the fraction of obscured AGN and its possible de- 
pendence on luminosity, by computing the relative frac- 
tion of IR emission in unobscured AGN spanning a large 
of luminosities. Using a sample of 206 optically-selected 
high luminosity AGN from the SDSS and 24 X-ray se- 
lected lower luminosity AGN from the GOODS and COS- 



MOS surveys with 0.8<z<1.2, nearly all having Spitzer 
detections at 24 /im, we found a decrease of a factor of 
~3 in the relative IR emission at rest-frame 12 fj,m from 
^boi=10 44 ergs s~ 4 to 10 47 5 ergs s _1 . In addition, we 
found a significant scatter of a factor of 2-3 at a given 
bolometric luminosity. Some of this scatter can likely be 
explained by a dependence of the 11-13 /zm fraction /12 
on viewing angle; some may also be due to an intrinsic 
scatter in / bsc at fixed bolometric luminosity due to a 
dependence on other parameters such as black hole mass, 
magnetic forces, etc. 

Using IR re-emission models in order to convert 
^Mips/iboi into a fraction of obscured AGN, we found 
that / bsc changes from ~ 90% at Lboi = 10 44 ergs s _1 to 
~ 30-40% at Lboi ^ 10 47 ergs s _1 . The derived depen- 
dence of the obscured fraction on luminosity is in good 
agreement with the direct, but possibly biased (by the 
omission of Compton-thick objects) observation of this 
fraction in X-ray surveys. We caution, however, that 
there may be systematic error in the precise values of 
/obsc inferred from iMiPs/^boi because our estimate of 
(/12) is likely to depend somewhat on the details of the 
density distribution within the torus. Nonetheless, it is 
interesting that in general the values of /obsc we infer are 
larger than what is seen in X-ray surveys. This contrast 
hints that there may be a significant population of heav- 
ily obscured, even Compton thick, AGN that are missed 
in X-ray observations but included in our samples. 
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TABLE 1 

Observed Properties of Sources in the SDSS Sample 



Name Redshift Optical Mag. (AB) Flux (10 -12 erg cm" 2 s" 1 ) log(Bol. Lurr 

u g r i z 24 /»m NUV FUV erg s" 1 

flux flux upper lower flux upper lower 

SDSS J110116. 41+572850.5 0.80 17.85 17.53 17.48 17.46 17.40 0.331 1.335 1.341 1.332 0.466 0.470 0.462 46.5 

SDSS J142810.31+353847.0 0.80 19.09 18.79 18.74 18.89 18.71 0.120 0.431 0.436 0.426 0.179 46.0 

SDSS J104755.02+120850.2 0.81 19.84 19.43 19.16 19.02 18.71 0.537 0.202 0.214 0.191 0.012 0.016 0.008 45.8 

SDSS J143345. 10+345939.9 0.81 20.36 19.64 19.10 18.84 18.49 0.485 0.089 0.091 0.086 0.055 45.5 

SDSS J160630.60+542007.5 0.82 19.03 18.77 18.73 18.83 18.61 0.413 0.455 0.459 0.452 0.185 46.1 

SDSS J235948. 53-103938. 5 0.83 18.07 17.87 17.79 17.80 17.71 0.439 1.209 1.226 1.190 0.421 0.436 0.407 46.5 



Note. — This table is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance 
regarding its form and content. 



TABLE 2 

Observed Properties of Sources in the GOODS Sample 
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TABLE 3 

Observed Properties of Sources in the COSMOS Sample 



Name Redshift Optical Mag. (AB) Flux (10~ 12 erg cm -2 s" 1 ) log(Bol. 
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